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"P^ . We give the nonperturbative phase diagram of the four-dimensional hot electroweak phase 

transition. The Monte-Carlo analysis is done on lattices with different lattice spacings (a). A 
p ■ systematic extrapolation o — > is done. Our results show that the finite temperature SU(2)- 

Higgs phase transition is of first order for Higgs-boson masses < 66.5 ± 1.4 GeV. At this 
endpoint the phase transition is of second order, whereas above it only a rapid cross-over can be 
seen. The full four-dimensional result agrees completely with that of the dimensional reduction 
^ ■ approximation. This fact is of particular importance, because it indicates that the fermionic 

sector of the Standard Model can be included perturbatively. We obtain that the Higgs-boson 
endpoint mass in the Standard Model is 72.4 ± 1.7 GeV. Taking into account the LEP Higgs- 
boson mass lower bound excludes any electroweak phase transition in the Standard Model. 
PACS Numbers: 11.10.Wx,11.15.Ha 



The observed baryon asymmetry is finally determined at the electroweak phase transition 
(EWPT) [|TJ. The understanding of this asymmetry needs a quantitative description of this phase 
transition. Unfortunately, the perturbative approach breaks down for the physically allowed Higgs- 
boson masses (e.g. m H > 70 GeV) 0. In order to understand this nonperturbative phenomenon a 
systematically controllable technique is used, namely lattice Monte-Carlo (MC) simulations. Since 
merely the bosonic sector is responsible for the bad perturbative features (due to infrared problems) 
the simulations are done without the inclusion of fermions. The first results dedicated to this ques- 
tions were obtained on four-dimensional lattices ||. Soon after, simulations of the reduced model 
in three-dimensions were initiated, as another approach [[|. This technique contains two steps. The 
first is a perturbative reduction of the original four-dimensional model to a three-dimensional one by 
integrating out the heavy degrees of freedom. The second step is the nonperturbative analysis of the 
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three-dimensional model on the lattice, which is less CPU-time consuming than the MC simulation 
in the four- dimensional model. The comparison of the results obtained by the two techniques is not 
only a useful cross-check on the perturbative reduction procedure but also a necessity. The reason for 
that is that the fermions, which behave as the heavy bosonic modes, must be included perturbatively, 
anyhow. 

In the recent years exhaustive studies have been carried out both in the four- dimensional || and 
in the three-dimensional |J sectors of the problem. These works determined several cosmologically 
important quantities such as the critical temperature (T c ), interface tension (a) and latent heat (Ae). 

Previous works show that the strength of the first order EWPT gets weaker as the mass of the 
Higgs-boson increases. Actually the line of the first order phase transitions, separating the symmetric 
and broken phases on the mg — T c plane has an endpoint, m# iC . There are several direct and indirect 
evidences for that. In four dimension at m# ~ 80 GeV the EWPT turned out to be extremely 
weak, even consistent with the no phase transition scenario on the 1.5-cx level [[7]. Three-dimensional 
results show that for tuh > 95 GeV no first order phase transition exists || and more specifically 
that the endpoint is rriH,c ~ 67 GeV 0. In this letter we present the analysis of the endpoint on 
four dimensional lattices. We study the thermodynamical limit of the first Lee- Yang zeros of the 
partition function |§. In order to get rid of the finite lattice spacing effects a careful extrapolation 
to the continuum limit is performed. The endpoint value of the SU(2)-Higgs model is perturbatively 
transformed to the full Standard Model (SM). 

We will study the four-dimensional SU(2)-Higgs lattice model on asymmetric lattices, i.e. lattices 
with different spacings in temporal (a t ) and spatial (a s ) directions. Equal lattice spacings are used 
in the three spatial directions (a, = a s , i = 1,2,3) and another one in the temporal direction 
(04 = a t ). The asymmetry of the lattice spacings is given by the asymmetry factor £ = a s /a t . The 
different lattice spacings can be ensured by different coupling strengths in the action for time-like 
and space-like directions. The action reads 

spM = 

A E f 1 - i Tr u v) + A E f 1 - \^ u pl 

sp ^ Z ' tp ^ Z 
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+ ^-Tr(^) + A -Tr(<p+<p x )-l 
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-K a Tr (wt+tPx* <Px) - ^Tr (<p+ +i U xA <f a 
fj,=i 

where U Xjfi denotes the SU(2) gauge link variable, U sp and U tp the path-ordered product of the four 
U x> n around a space-space or space-time plaquette, respectively; <p x stands for the Higgs field. It is 
useful to introduce the hopping parameter and f3 2 = (3 s (3 t . The anisotropies 7? = f3 t / (3 S and 

7k = K t/ K s are functions of the asymmetry £. These functions have been determined perturbatively 



10] and non-perturbatively [11] demanding the restoration of the rotational symmetry in different 
channels. In this paper we use the asymmetry parameter £ = 4.052, which gives r y K = 4 and 
7/3 = 3.919. Details of the simulation techniques can be found in p. 

We have performed our simulations on finer and finer lattices, moving along the lines of constant 
physics (LCP). In our case there are three bare parameters (k, (3, A). The bare parameters are chosen 
in a way that the zero temperature renormalized gauge coupling is held constant and the mass 
ratio for the Higgs- and W-bosons Rhw = ran/raw corresponds to the Higgs mass at the endpoint 
of first order phase transitions: Rhw,c- These two conditions determine a LCP as a one-dimensional 
subspace in the original space of bare parameters. The position on the LCP gives the lattice spacing 
a. As the lattice spacing decreases Rhw,c R-HW,amt.- A schematic illustration is shown in Fig. 1. 
The LPC (solid line) defined by the endpoint represents the above idea. The short dashed lines give 
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Figure 1: Schematic view of the phase diagram. The solid line represents the LCP defined by the 
endpoint condition. The numbers on the line correspond to the temporal extension for which the 
endpoint is realized (the dashed lines show their projection to the k - A plane). The dotted lines 
running into these numbered points correspond to first order phase transitions for g 2 R = const, but 
different Rhw~s- A LCP defined by a constant Rhw value is shown by the long dashed line. 



the projections to the A — k plane. The increasing numbers on the LCP show the temporal extensions 
of the lattice, thus corresponding to smaller and smaller lattice spacings. The dotted lines represent 
phase transition points of theories with fixed renormalized g 2 and L t but different Rhw values. Along 
the dotted lines one can observe first order phase transitions upto the LCP defined by the endpoint 
condition. Note, however, that this endpoint LCP is not the same as the LCP defined by the constant 
Rhw — Rnw,cont. value (long dashed line). They merge for decreasing lattice spacings, but at larger 
a the difference is the result of the 'poor realization' of Wilson's RG transformations with only three 
terms and parameters in the action. It is worth mentioning that the SU(2)-Higgs model is trivial 
for small gauge couplings, therefore, the a —>■ limit can not be performed. Even the points on the 
endpoint LCP do not define continuum theories. The second order phase transitions on it merely 
reflect a finite temperature phenomenon, the corresponding zero temperature SU(2)-Higgs theory is 
still trivial. 

Since our theory is a bosonic one we assumed that the finite size corrections are quadratic in the 
lattice spacings; therefore an a 2 fit has been performed for Rhw,c in order to determine its continuum 
value. 

The technical implementation of the above LCP idea has been done as follows. By fixing (3 = 8.0 
in the simulations, we have observed that is essentially constant within our errors. For the small 
differences in g^ we have performed perturbative corrections. We have carried out T ^ simulations 
on L t = 2, 3, 4, 5 lattices (for the finite temperature case one uses L t <C L x , L y , L z ), and tuned k to 
the transition point. This condition fixes the lattice spacings: at = a s /C, = l/(T c L t ) in terms of the 
transition temperature T c in physical units. The third parameter A, finally specifying the physical 
Higgs mass in lattice units, has been chosen in a way that the transition corresponds to the endpoint 
of the first order phase transition subspace. 

In this paper V = L t ■ L\ type four- dimensional lattices are used. For each L t we had 8 different 
lattices, each of them had approximately twice as large lattice-volume as the previous one. The 
smallest lattice was V = 2-5 3 and the largest one was V = 5-50 3 . We collected quite a large statistics 
and the Ferrenberg-Swendsen reweighting |Tj| was used to obtain information in the vicinity of a 
simulation point. 
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0.000178 


0.107733 


0.934(10) 


0.569(4) 


3 


0.000178 


0.106988 


0.913(12) 


0.575(3) 


4 


0.000178 


0.106620 


0.905(8) 


0.585(5) 


5 


0.000178 


0.1064974 


0.867(36) 


0.566(30) 



Table 1: Summary of simulation parameters and results on Rhw and g% at T = 0. 



The determination of the endpoint of the finite temperature EWPT, thus a characteristic feature 



of the phase diagram, is done by the use of the Lee- Yang zeros of the partition function Z |[3] . Near 
the first order phase transition point the partition function reads 



Z = Z s + Z b oc exp(-Vf s ) + exp(-Vf b 



(2) 



where the indices s(b) refer to the symmetric (broken) phase and / stands for the free-energy densities. 
Near the phase transition point we also have 



fb = fs + a(n - k c ) , 
since the free-energy density is continuous. It follows that 

Z = 2exp[-V(f s + f b )/2] cosh[-Va( K - k c )} 
which shows that for complex k Z = vanishes at 

Im(re) = 7T ■ (n - 1/2) /(Va) 



(3) 



(4) 



(5) 



for integer n. In case a first order phase transition is present, these Lee- Yang zeros move to the 
real axis as the volume goes to infinity. In case a phase transition is absent the Lee- Yang zeros stay 
away from the real k axis. Thus the way the Lee- Yang zeros move in this limit is a good indicator 
for the presence or absence of a first order phase transition ]n| Denoting k$ the lowest zero of Z, 
i.e. the position of the zero closest zero to the real axis, one expects in the vicinity of the endpoint 
the scaling law Ihi(k ) = ci(L t , \)V U + c 2 (L t , A). In order to pin down the endpoint we are looking 
for a A value for which c 2 vanish. In practice we analytically continue Z to complex values of k 
by reweighting the available data. Also small changes in A have been done by reweighting. As an 
example, the dependence of c 2 on A for L t = 3 is shown in fig. 2. To determine the critical value of A 
i.e. the largest value, where c 2 = 0, we have performed fits linear in A to the nonnegative c 2 values. 

Having determined the endpoint A crit .(L f ) for each L t we calculate the T = quantities (Rhw, 9r) 
on V = (32L t ) ■ (8L t ) ■ (QL t ) 3 lattices, where 32L t belongs to the temporal extension, and extrapolate to 
the continuum limit. All the T = simulations were performed at A = 0.000178 and an extrapolation 
to the X cr it.(Lt) has been made. The parameters and results of the simulations are collected in 
table 1, while table 2 shows the Rhw values extrapolated to the A cr j t .(L t ). Having established the 
correspondence between X crit (L t ) and Rhw, the L t dependence of the critical Rhw is easily obtained. 
Fig. 3 shows the dependence of the endpoint Rhw values on 1/Lj. A linear extrapolation in 1/Lf 
yields the infinite volume (i.e. continuum limit) value of the endpoint Rhw- We obtain 66.5 ± 1.4 
GeV, which is our final result. 

Comparing out result to those of the 3d analyses one observes complete agreement. Since 
the error bars on the endpoint determinations are on the few percent level, the uncertainty of the 
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2 


0.0001773(14) 


0.1077292(2) 


0.932(10) 


3 


0.0001664(27) 


0.1069581(2) 


0.883(12) 


4 


0.0001590(44) 


0.1066316(3) 


0.856(8) 


5 


0.0001664(20) 


0.1064948(6) 


0.838(36) 



Table 2: Critical A corresponding to the endpoint of phase transition as function of L t and the 
corresponding value of Rhw- 




X * 10 4 

Figure 2: Dependence of c 2 on A for L t = 3. 



dimensional reduction procedure is also in this range. This indicates that the analogous perturbative 
inclusion of the fermionic sector results also in few percent error on M#. 

Based on our published data || |llj and the results of this paper we are now able to draw the 
precise phase diagram of the SU(2)-Higgs model in the {T c /m H — Rhw) plane. This is shown in fig. 
4. The continuous line - representing the phase-boundary - is a quadratic fit to the data points. 

Finally, we determine what is the endpoint value in the full SM. Our nonperturbative analysis 
shows that the perturbative integration of the heavy modes is correct within our error bars. Therefore 
we use perturbation theory JL4| to transform the SU(2)-Higgs model endpoint value to the full SM. 
We obtain 72.4 ± 1.7 GeV, where the error includes the measured error of RHW,cont., g% and the 
estimated uncertainty [|15j due to the different definitions of the gauge couplings between this paper 



and fTM. The dominant error comes from the uncertainty on the position of the endpoint. 

In conclusion, we have determined the endpoint of hot EWPT with the technique of Lee- Yang 
zeros from simulations in four- dimensional SU(2)-Higgs model. The phase diagram has been also 
presented. The phase transition is first order for Higgs masses less than 66.5±1.4 GeV, while for larger 
Higgs masses only a rapid cross-over is expected. One of the most important results of the present 
letter is that integrating out the heavy modes perturbatively is precise as shown by a comparison to 
our nonperturbative results. Thus the above 66.5 ±1.4 GeV value can be perturbatively transformed 
to the full SM. We obtain 72.4 ±1.7 GeV for the endpoint Higgs mass. As pointed out above the 
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Figure 3: Dependence of Rhw,c, i-e. Rhw corresponding to the endpoint of first order phase 
transitions on 1/Lf and extrapolation to the infinite volume limit. 



perturbative inclusion of the fermionic sector of the SM is also correct to a few percent error level. 



The present experimental lower limit of the SM Higgs-boson mass is 89.8 GeV | lq| . Taking 
into account all errors (in particular those coming from integrating out the heavy fermionic modes), 
our endpoint value excludes the possibility of any EWPT in the SM. This also means that the SM 
baryogenesis in the early Universe is ruled out. 

More details of this investigation will be published in a forthcoming publication ]17 . 
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Figure 4: Phase diagram of the SU(2)-Higgs model in the (T c /mH — Rhw ) plane. The continuous 
line - representing the phase-boundary - is a quadratic fit to the data points. 
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